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t^J- ■ Abstract 

t— i ; 

In recent time, by working in a plane with the metric associated with wave equation (the Spe- 
cial Relativity non-definite quadratic form) , a complete formalization of space-time trigonometry 
Q_r and a Cauchy-like integral formula have been obtained. 

^ fi-i In this paper the concept that the solution of a mathematical problem is simplified by using a 

"mathematics" with the symmetries of the problem, actuates us for studying the wave equation 
(in particular "the initial values problem" ) in a plane where the geometry is the one "generated" 
i— — i. by the wave equation itself. 

In this way, following a classical approach, we point out the well known differences with 
■ respect to Laplace equation notwithstanding their formal equivalence (partial differential equa- 

^sO . tions of second order with constant coefficients) and also show that the same conditions stated 

for Laplace equation allow us to find a new solution. In particular taking as "initial data" 
£T) " for the wave equation an arbitrary function given on an arm of an equilateral hyperbola, a 

. "Poisson-like" integral formula holds. 

o: 

Keywords: Hyperbolic geometry. Wave equation. Boundary value problem. 
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X 1 

^- 1 Introduction 

In a recent paper [1] and books [2], [3] it has been shown how a complete formalization of Minkowski's 
space-time geometry and trigonometry has been obtained by means of hyperbolic numbers that, from 
an algebraic point of view, are the simplest extension of complex numbers and are defined as 

{z = x + hy; h 2 = 1 ; x, y G R; h g R}. 

These results have been obtained working in a Cartesian plane with the non-definite metric corre- 
sponding to the modulus of hyperbolic numbers (square distance: d 2 = x 2 — y 2 ), so as the Euclidean 
distance corresponds to modulus of complex numbers [3]. Here we call hyperbolic the plane with this 
metric. 
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By means of this approach, in [2] the bases for studying the functions of a hyperbolic variable 
have been set and, for these functions, a Cauchy-like integral formula has been stated [4]. 

Since these results refers to functions satisfying the two-dimensional wave equation here we begin 
to investigate if for studying this equation, in particular with regard to the initial value problem, it 
is appropriate to work in the hyperbolic plane. 

Practically we came back from the results of Special Relativity in the following meaning: 
Invariance of wave equations =>- Lorentz transformations =>- Minkowski (hyperbolic) geometry =>- wave 
equation in hyperbolic plane. 

In this approach, save for the recalled novelty, we follow Riemann who obtained his integral formula 
[5, p. 450] as a precursory of special relativity in the meaning that for studying the initial value 
problem for the wave equation, he considered as equivalent space and time and represented them as 
coordinates in a plane. 

In the appendix we "translate" in the hyperbolic plane some properties that hold in Euclidean 
geometry. 

2 Initial Data Problem for Wave Equation 

Following the Cauchy theory [5], the solution of a partial differential equation (PDE) of degree N 
can be obtained by a series development around a point in which the values of the function and its 
partial derivative, up to degree N — 1, are given. 

As an exception to this approach, the solution of the second degree partial differential Laplace 
equation, is determined by just the values of one arbitrary function given on the frontier of a domain 
with "appropriate regularity conditions" [5, Chap. IV]. This problem is known as Dirichelet's problem 



As a difference from Laplace equation for which the solution is determined inside a closed domain, 
for the wave equation the initial data are given on an open, appropriate curve [5] and the solution is 
determined in the two opposite sides of the given curve, in particular in the domain determined [5, 
p. 450] by the parallel to axes bisectors (characteristic lines) from the extreme points of the curve. 
In this paper we follow the Euclidean approach to Laplace equation and translate it to a "hyperbolic 
approach" to wave equation. We see that, studying the problem in this way, a new situation arises. 
In particular, by means of a "Poisson-like" integral, we obtain the sum of the solutions in two points 
that we define as symmetric by extending to hyperbolic plane the well known symmetry with respect 
to a circle in Euclidean plane [2, Sect. 7.5]. 

We begin by translating the basic mathematics into hyperbolic geometry. 

2.1 Integral Identity in Hyperbolic Plane 

For studying the initial data problem for partial differential equations, the Green's formulas are 
usually applied [5]. They represent a particular application of Gauss formula 



that transforms the integral on a domain into an integral on its frontier. 

Let us now apply this identity to wave equation and look for its application to the "initial value 
problem" following the Green's approach. 



[5, Chap. IV §2]. 




(1) 
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Let us consider two arbitrary functions u, v continuous, with the derivatives that appear in the 
formulas, in a domain D up to its contour T [5, Chap. IV]. Let us set 

du ^ r du 

= "a? Y = (2) 

and introduce the differential parameters [2, Chap 8] in the flat hyperbolic plane, i.e., with pseudo- 
Euclidean metric 

A ^ U ^ U ( b ) Af dv du dv du 

2 dx 2 dy 2 ' ' dx dx dy dy 

Equation (1) becomes 

J j v A 2 u dx dy + J j A(v , u) dx dy = (j) v ^— d y + — d x^j . (4) 

Now we go on as it is usually done for studying the "initial values problem" for Laplace equation. 
By subtracting from Eq. (4) the equation obtained by changing v ■<->■ u, we obtain a Green identity 
for the differential operators (3 a) 



J J (v A 2 u — u A 2 v ) dx dy = <j> 



du du \ f dv dv 
— dy + — dx - u — dy + — dx 
dx dy J \dx dy 



(5) 



In order to calculate the line integral along the curve T, on the right hand side of Eq. (5), a local 
reference frame is used. This frame has its origin in the point that moves along the curve and an 
axis (r) tangent to the curve, oriented according with the integration direction. 
For the construction of the other local axis n, the hyperbolic geometry is applied to the plane x, y. 
So we take the n axis in the direction of the hyperbolic normal to r axis and oriented so that the 
frame n, r is congruent with the orientation of x, y frame. According with the topology of hyperbolic 
plane [1], different kinds of pairs of unity vectors originate. 

A detailed treatment of this subject, based on [1], is developed in App. A and summarized in the 
caption of Fig. 3 where the hyperbola with \dy/dx\ > 1, that is the curve employed in this paper, is 
considered. 

Now we write the argument in brackets of the integrals in the right hand side of Eq. (5) as a function 
of the local coordinates n, r. 
Therefore, by setting 

• du/dn the derivative of u in the direction of n; 
taking into account that 

• in the line integration it results dn — 0, 

• the relations between the derivatives with respect to the orthogonal directions of the tangent 
and the "hyperbolic normal" to a curve are (Eq. (95)) 

dx dy dx dy 
dr dn ' dn dr' 
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we have 



du du , f du dy dudx\ f du dx dudy\ du . . 

d v + ^r dx = «r 17- + 77717- ^= 777 777: + 777777/ dr = i^ dr - ( 7 ) 



cte <% \9x <9r dy dr ) \dx dn dy dn J dn 

With these definitions and expressions Eq. (5) becomes 



/ f (v A 2 u — u A 2 v) dx dy — (f ( v ^ — u tt-) dr. 
JdJ ~ ' Jr\ dn dn) 



(8) 



Therefore, by working in hyperbolic plane, i.e., with the geometry related with wave equation [2], 
the same integral identity (8) that holds by applying Euclidean geometry to x, y plane in the study of 
Laplace equation [5, p. 252, Eq. (26)], has been obtained. 

Now we go on as it is usually done for studying the "initial values problem" for Laplace equation. 
In particular if u, v satisfy the wave equation, the left-hand side of Eq. (8) is zero and we have 

dr = a (9) 

dn dn I 



By setting in Eq. (8) v — 1 we obtain 



D 



du 

A 2 u dxdy = <p — dr, (10) 



and if the function u satisfy in D the wave equation and holds the appropriate regularity conditions, 
that for Laplace equation are: to be continuous with the partial derivatives in the domain and satisfy 
the Holder conditions on the contour [6], [7, p. 50], we have: 

The integral on the contour (T) of the derivative of u with respect to the normal is zero: 

du , , . 

s <fr = 0. (11) 



2.2 Characteristic Domains in Hyperbolic Plane 

In the application of the relations of the previous section to the "initial value problem" for Laplace 
equation and in complex analysis for demonstrating the Cauchy integral formula, the singularity at 
a point Q = (x, y) is excluded from the domain of integration by means of a circle, centered in Q, 
with radius r — > 0. In Euclidean plane this circle represents the locus of points at the same distance 
from Q and allows one to obtain useful simplifications and relevant results. 

In hyperbolic plane the locus of points at the same distance from a given point is the equilateral 
hyperbola, then in [4] the circle with radius r — > has been replaced, by an equilateral hyperbola 
with semi-diameter p — > 0. This means that the hyperbola becomes the parallel to axes bisectors 
from the point Q. 

This domain can be recognized as the one considered by Riemann for his approach to the initial 
value problem for wave equation [5, p. 450]. 

Moreover for Laplace equation the "initial data" are given on a closed domain (a topological 
transformation of circles), for the wave equation, they are given on a curve that can be considered 
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as a topological transformation of an equilateral hyperbola, i.e., of a curve with tangent lines of a 
given kind [5]. 

Now we observe that for the "initial data" for wave equation given on a curve, we can consider, 
from a mathematical point of view, both the points at the left or at the right of the curve. This 
possibility generate the difference [5] between the Laplace and wave equations with respect to the 
initial data problem. In particular, as we see in the following of this paper, for one arbitrary function 
given on an arm of equilateral hyperbola, we do not have an unique determination of a function 
satisfying the wave equation and assuming on the line the given values, as it happens for Laplace 
equation, but from these values we can determine the sum of the values of the function in two points 
that, extending to hyperbolic plane the Euclidean symmetry about the circle, can be defined as 
symmetric with respect to hyperbola (Sect. 3.1). 

Taking into account the analogies and these differences, we now "translate" to wave equation, 
studied in the hyperbolic geometry, the classical results obtained for Laplace equation, studied in 
the Euclidean geometry, in the internal points of a circle [5]. 



3 Initial Data on an Arm of Equilateral Hyperbola 

Let be given, as initial data, the values of one arbitrary function u (t) on the arc, between the points 
A and B, of the right arm of an equilateral hyperbola 7 with center in the axes origin O and semi- 
diameter p (Fig. 1). 

We look for a function u ( x, y) satisfying the wave equation and assuming on 7 the given values. 
By calling Q = (£, 77) the point in which we look for u(£, rj), the parallels to axes bisectors from Q 
cross 7 in the points Pi , P 2 . These points have to be internal to hyperbola arc AB. 
In this way the domain of dependence, determined by the point Q and the hyperbola 7, is defined 
[5, p. 438]. 

We draw from P 1 and P2 the parallel to axes bisectors in the opposite direction with respect to 
Q and call Q* their intersection point. 

We define Q and Q* as symmetric with respect to hyperbola 7, by extending to hyperbolas in the 
hyperbolic plane, the Euclidean symmetry with respect to a circle. Actually these points, as it is 
shown in [4], have equivalent properties: 

1. they are on the same straight line through the center of hyperbola 7; 

2. the product of their hyperbolic distances from the center is equal to p 2 (squared semi-diameter). 
Let us set Q* = (£*, r]*) and call 



it results 



OQ = ^e^f = q, (12) 

C = ^Z, V* = ^V- (13) 

The other two hyperbolas X and X* with centers in Q and Q*, represented in Fig. 1, are taken 
so that they intersect each other on the hyperbola 7 in the points P{ and P 2 \ We call their semi- 
diameters p and p*, respectively. 

For the other elements of Fig. 1, we call a the hyperbolic angle between x axis and the straight line 
OQ and 0, 9, 9* the angles describing the hyperbolas, measured with respect to the straight line OQ. 
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The equations of the three hyperbolas, in hyperbolic polar form, are given by 

7 — » x = p cosh(0 + a) ; y = p sinh(0 + a) (14) 
X a; = £ + p cosh(# + a) ; y — r] + p sinh(# + a) , (15) 

X* x = C - p* cosh(0* + a) ; y = r]* - p* sinh(0* + a) . (16) 

Moreover we consider a point P = (x, y) and the distances r = QP and r* = Q*P given by 

r = v /( x _^)2_ (y _ r? )2. r* = v/(x-e) 2 "(y-^) 2 , (17) 

where £* and ?7* are given by Eqs. (13). 

In Fig. 2 we report the elements of Fig. 1 rotated by a hyperbolic angle —a. 
In this way a symmetric representation with respect to x axis is obtained. This representation allows 
a better insight in the geometric properties and an easier way for definitions and calculations. 

After this rotation, we have 

• the hyperbola 7 remains in the same position; 

• all the hyperbolic distances and hyperbolic angles between corresponding lines are preserved; 

• the hyperbolas X and X* become symmetric with respect to x axis. 

Thus, by calling 0*, 9 l , and 9 l * the absolute values of the limit hyperbolic angles for a point moving 
along the aforesaid hyperbolas from P{ to P\, the ranges for the hyperbolic angles in Eqs. (14-16), 
are 

for 7 : < < , for X : -d l < < +9 i , for X* : +9 { * > 9* > -9" . (18) 



3.1 Extension of Apollonius Circle Theorem to Hyperbolic Plane 

We show that a similar theorem to the Apollonius theorem in Euclidean plane about the circle, holds 
for equilateral hyperbolas in hyperbolic plane. 
Referring to letters and symbols used in Fig. 2, we have 

Theorem - Given a straight line from the center of a hyperbola with semi- diameter p and two points 
Q and Q* on this straight line, so that 



35.35^22 = ^, (19) 



for all the points P on the hyperbola, we have 



2L = 99.. ( 2 o) 

Q*P p v 1 

This theorem can be demonstrated in a "Euclidean way" thanks to the analytical formalization 

of Hyperbolic geometry [1]. Here we use the analytical approach that has allowed the recalled 
formalization. 

Proof - By using the definitions of the hyperbolic distances r and r* given by Eq. (17) and the 
definition of q given by Eq. (12), Eq. (20) becomes 



r _ V( x ~ 2 - (y ~ v) 2 _ 1 ^ 



r * v/(x - e) 2 - (y ~ V*) 2 P 
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By squaring this equation, it results 

( r Y- {x 2 -y 2 ) + (j 2 -n 2 )-2(jx-ny) _ q 2 
\r* J (x 2 — y 2 ) + (£* 2 — n* 2 ) — 2(£* x — n* y) p 2 

By means of Eqs. (12) and (13), we see that Eq. (22) is verified if P = (x, y) belongs to 7, that is 
x 2 — y 2 = p 2 . Actually, in this case, we have 



2 



r\ 2 _ p 2 + q 2 — 2(1; x — r/y) q 
r *J - p 4 2 p 2 
P + — -2— itx-rjy) 



□ (23) 



For convenience we set 



and Eq. (21) becomes 



1 A p (24) 



p 



(4) = A p- (25) 



P on 7 

From this theorem, by applying Eq. (25) in the particular cases P = P[ or P = P 2 \ it follows 



p _ QP{ _ QP\ _ ( r 



P PIQ* P\Q* Vr*/p „ 7 



= A P . (26) 



3.2 Application of Integral Formulas 

Let us apply Eq. (9) to the following domains represented in Fig. 1 

1. between the hyperbolas X and 7; 

2. between the hyperbolas 7 and X*. 

Let u ( x, y ) be a function that satisfies the wave equation in these domains. 
In the application of Eq. (9) to the domain 1 we set 

v(x, y) = lnr , (27) 

in the application to the domain 2 we set 

v(x,y) = \nr% (28) 

where r, r* are given by Eq. (17). 

It can be checked at once that these functions v(x,y) satisfy the wave equation. 
From Eq. (9) we have: 
for domain 1 



[ p * ( Ou dlnr\ ^ [ p i ( Ou dlnr\ , n , nft , 
/ lnr— -u— — ) dr+ / lnr — -u— — Wr = (29) 
ip 1 i ( 7 ) V dn dn J Jp\(x) \ on dn J 



for domain 2 



p 2 i ( 7 ) V ^ n dn J jp{(t*) V ^ n ^ n 



7 




Figure 1: The integration domains for wave equation in hyperbolic geometry. 

In this figure geometric elements are represented for application of integral formulas to calculate a 
function satisfying the wave equation, in the symmetric points Q (£, rj) and Q* (£*, n* ), with initial 
data given by an arbitrary function defined on the arc AB of an equilateral hyperbola 7. 
Pi and P 2 are the extreme points of the domain of dependence ( Sect. 3). In particular the following 
elements are reported 

• The equilateral hyperbolas 7, X, X*, with semi-diameters p, p, p* and their intersection points PI 
and P*. 

• The hyperbolic angular variables 0, 9, 9*, measured with respect to the straight line OQ*, set at a 
hyperbolic angle a with respect to x axis. 

• The hyperbolic distances r and r* from a point P of the hyperbola 7, to points Q and Q*, 
respectively. 

Since the hyperbolas, in hyperbolic geometry, represent the locus of points at the same distance from 
a given point, they correspond to circles of Euclidean geometry, used for the same problems about 
Laplace equation studied in Euclidean plane. 

In particular, the hyperbolas X and X*, for which, in the final step of the procedure, we do the limit 
p, p* — > 0, correspond to the infinitesimal circles around the singularities. In this limit they become 
the parallel to axes bisectors from the points Q and Q*, respectively. 
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Figure 2: Domain of dependence in symmetric position. 

The coordinates of the points of Fig. 1 are transformed by means of a hyperbolic rotation of an angle 
— a in order to set Q on the x axis. In this way a symmetric representation with respect to x axis is 
obtained. This representation allows a better insight in the geometric properties and an easier way 
for definitions and calculations. 
After this rotation, we have 

• the hyperbola 7 remains in the same position; 

• all the hyperbolic distances and hyperbolic angles between corresponding lines are preserved; 

• the hyperbola X and X* become symmetric with respect to x axis. 

• the extreme points A and B of the arc on which the initial data u (r), are given, do not change. 
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Let us add Eqs. (29) and (30) and group together the integrals on 7 and the integrals on X and X* . 



As far as the integrals on 7 are concerned, for Eq. (25), it results In r* 
we have 



In r — In A„ . Therefore 



J p\ 



2 1 9U A 

m r — — a r 

(7) dn 



f 



din r 
u — — dr + 

on 



■ lnAp 7h; dr ' 

j p\ ( 7 ) ° n 



r 

J Pi (■ 



1 9U A 

in r — a r 

on 



2(7) 

Pi dlnr* , 
u — dr . 

Pi (7) 



dn 



(31) 



The first and third integrals are equal, but have opposite integration direction then their sum is zero. 
So, by collecting the second and the fifth integrals, Eq. (31) becomes 



n 1 a 9u a r 

/ In A p — d t - / 

J P i ( 7 ) on J P i ( 7 ) 



u 



d In (r/r*) 
dn 



dr. 



(32) 



Now let us consider the integrals on X and X*. They have to be calculated for r = p and r* = p*, 
respectively. From the relation (26), we have In p* = In p — In A p , therefore we have 



JPl(X) on J Pi(x) \ dn ) {r=p) J F 

pi 



du , 
In p -^r- dr + 



+ / lnA p — dr - / 

Jpi(i*) an J p\ 



u 



pi (i*) 

d In r" 

dn 



dn 



dr . 



(33) 

The sum of the first and third terms is an integral on a closed cycle and, because lnp is constant 
and u satisfies the wave equation, it follows, from Eq. (11), that it is zero. For the same reason the 
integral on the closed cycle given by the sum between the first term of Eq. (32) and the fourth one 
of Eq. (33) is equal to zero. 

After these reductions, the contribution of the derivative du/dn disappears and from the sum of Eqs. 
(29) and (30) remains 



r " 

j pirn 



din r 

dn 



dr ■ 1 11 

(r=p) J Pi (X*) 



d In r* 

dn 



dr 



(r*=p*) 



U 



P{ (7) 



d In (r/r*) 
dn 



dr. 



(34) 



3.3 Introduction and Properties of C Function 

Let us now introduce a function C given by 



In 



E (4) 

a \r* / 



By using this function, the right hand side of Eq. (34) may be written 



u 



PI (7) 



d In (r/r* 
dn 



dr = 



pi 



pi (7) 



u — dr . 

dn 



(35) 



(36) 



The C function has the following properties 

• thanks to relation (25), it is zero on 7, 

• it satisfy the wave equation. 

Therefore it has the same properties of Green function introduced for Laplace equation [5, Chap. 4]. 
Here we see that, if the "initial data" are given on an arm of equilateral hyperbola, the C function, 
as the Green function, allows us to obtain a Poisson-like integral. 
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4 Poisson-like Integral 

The Poisson integral formula originates for solving in a circle (and in a sphere), the "initial data" 
problem for two (and three) dimensional Laplace equation and states a different way for tackle the 
problem for elliptic partial differential equations (PDE), with respect to Cauchy problem, about 
initial data (Sect. 2). It is known as "Dirichelet problem" [8, p. 345]. For two-dimensional Laplace 
equation it may be obtained in many ways [5] - [8]. Here, studying the problem in hyperbolic plane, 
we obtain an analogous integral formula for "initial data" given on the right arm of equilateral 
hyperbola 7. We call it Poisson-like integral, for the wave equation. 



4.1 Poisson-like Kernel 



With reference to Eq. (34) let us consider the right hand side, in which the values of u are given. 
The first step is to calculate the kernel dC/dn of the integral. From Eq. (35) it results 



dC _ 1 dr 
dn r dn 



1 dr* 



r* dn 

Let us carry out the derivatives by considering r and r*, given by Eq. (17), as functions of x, y 

dr' 



dr dr dx dr dy 
dn dx dn dy dn' 



dr* dx dr* dy 
+ 



dn dx dn dy dn ' 
Moreover, because the point P = (x, y) is on 7, it results 



(37) 



(38) 



222 
x — y = p 



and, from Eqs. (25) and (97), the terms of the right hand side of Eq. (37) become 

1 



1 dr 

r dn 



- (y-v)- 
p p 



pr 



2 [P 2 - (£x-riy)] 



1 dr* 

r* dn 



( p 2 c ^ x i p 2 ^y 

q z p q l p 



pr A 



1 - \ {ix-ny) 



By substituting these expressions in relation (37), after reduction, we obtain 



dC_ 

dn 



on 7 



p 2 — q 2 
pr 2 



(39) 
(40) 

(41) 



Now we see that, by means of this expression, the integral in the right hand side of Eq. (34) gives a 
Poisson-like kernel. 

Proof - By referring to Fig. 1 let us calculate the hyperbolic distance r = QP. By applying the 

A 

hyperbolic Carnot's theorem [1] to the triangle QOP, where OP = p, OQ = q and is the hyperbolic 
angle QOP 

QP 2 = r 2 = p 2 + q 2 - 2pq cosh 0, (42) 
Moreover, on hyperbola 7 we have 

dr — pd(f) . 



(43) 
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Thus by means of Eqs. (41), (42) and (43), the right hand side of Eq. (34) becomes 



/ u{<P + a) f q —d<t>. □ (44) 

J-4fi{i) V + Q - 2pqcosh(p 

Referring to Fig. 1, we note that as points P{ and P 2 * go toward the extreme points Pi and P 2 , the 
sides QP[ and QP\ become parallel to axes bisectors and r — > 0. Therefore, in this limit, the integral 
diverges. 



4.2 Limit of Poisson-like Integral on 7 



By referring to Fig. 1, we can note that as P — > Pi and P 2 , more than r — > 0, we have 9 — > 00. This 
angle is linked with the integration variable 0. 

Here we see that the limits of the integral (44) is of the same order as 6 l — > 00. This fact allows us 
to obtain for r — > 0, a finite result. 

Proof - Let us divide the left and right sides of Eq. (34) by 9 l and let us begin by calculating the 
limit of the right-hand side 



lim 

6* -too 



2Q i 



+4> 1 



u{(p + a) — 
: ( 7) v + q 2 - 2 p q cosh <p 



(45) 



This limit is an indeterminate form 00/00. 

In order to calculate this limit, let us express the angle 6 l as function of <p\ that is the limit value of 
the integral. In particular for the point P\ = {p cosh0* ,p sinh0*) of Fig. 2, we have 



tanh^ = 



p sinh0j 



p cosh 4> l — q 
and, substituting in Eq. (45), we obtain 

( r+4> i 



lim < 



6 l =tanh 



p 2 — q 2 



-1 



p sinh ( 



p cosh 4> l — q 



(46) 



¥ (7) 



u ( (f) + a) — 

p 2 + q 2 — 2 p q cosh <p 



2 tanh -1 



p sinh i 



p cosh 0* — q 



> . 



(47) 



Let us apply the rule of L'Hospital by substituting to numerator and denominator their derivatives 
with respect to 0\ In this way the integral in the numerator is eliminated and we have 



u ( (j) 1 + a ) 



p 2 — q 2 



lim 



p 2 + q 2 — 2 p q cosh < 



+ u(-(f) i + a) 



p 2 — q 2 



p 2 + q 2 — 2 p q cosh(— c, 



p — pq cosh < 



(48) 



p 2 + q 2 — 2 p q cosh <p l 



By calculating cosh <p 2 , from the coordinates of extreme P 2 in Fig. 2, it results cosh0 2 = (p 2 + 
q 2 )/{2pq) and we have 

lim [«(<// + «) + w(-<^ + a)] =w(0 2 + a) + w(-0 2 + a). = u(P 1 )+u(P 2 ) □ (49) 
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Actually we can say that the limit of the ratio between the Poisson kernel and 2 9 l acts as a 
hyperbolic delta function. In fact the final result in Eq. (49) is that the integral disappears and just 
the sum of the values of the integrand calculated in the points Pi and P2 remains. These points are 
the ones connected, by the parallel to axes bisectors, with the points in which we are looking for the 
field. 



4.3 Limits of the Integrals on X and X* 

Now let us consider the two terms of the left hand side of Eq. (34) and calculate the same limits of 
the right hand side. 

Let us express the two integrals as functions of the hyperbolic angular variables 6 and 6* , respectively, 
and consider the local axis n in the points of X. It is directed like r and, referring to the right arm 
in Fig. 3, we observe that, because in Eq. (34) the integration direction is upward, n is oriented in 
the r increasing direction and dr/dn = 1. We have 

d In r ( d In r dr \ 1 

(50) 



dn \ dr dn J r=p p 
Moreover, in the points of X it results dr = p d6. 

From these positions and by recalling Eq. (18) that gives the ranges of the hyperbolic angles, the 
first term of the left hand side of Eq. (34) becomes 

' /<91nr\ f +e% 

u I — — — dr = u(p,9 + a)d9. (51) 



p\(t) \ dn /(r=p) j-0i(x) 
Analogous considerations allow us to transform the integral on X*. 

Let us consider the local axis n in the points of X*. It is directed like r* and, referring to the left 
arm in Fig. 3, we observe that, because the integration direction is downward, n is oriented in the 
r* increasing direction, then dr* /dn = 1 and we have 

dlnr* ( d In r* dr* \ 1 



dn \ dr* dn J r * =p * p* 
Moreover, in the points of X* it results dr = p* d6* . 

From these positions and by recalling Eq. (18), the second term of the left hand side of Eq. (34) 
becomes 

/ u — dr= u(p*,6* + a)d6* . (53) 

The values of 9 l and 9 l * are calculated from the intersection points of X and X* with 7, for a = 
(Fig. 2). By taking into account Eqs. (12) and (19) and setting 



q* = OQ* = V -, (54) 



it results 



cosh „< = P 2 -^-» 2 , sinhe . = vggWEgEgj , ,55) 

2qp 2qp ' 

coshe ^=l^Cl±fl, sinh ,,, = yTgHgF^SZEBg). (56) 

2q*p* 2q*p* v ' 
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Let us demonstrate, referring to the particular case of Fig. 2 (a = 0), that 

1 



lim 



lim 



29 



0i (X) 



u (p, 9) d9 



1 r ti 

V J -$**(!*) 



)d9* 



= u(Q) 
= u(Q*) 



(57) 
(58) 



Proof. Let us start from Eq. (57) and write the function u ( p, 9 ) on the points of hyperbola X 
by means of the Taylor's formula of order 1 



u(p, 9) = u(q + p cosh 6 1 ; p sinh 0) = u ( Q ) 



+ 



+ 2 



9m \ 



) p cosh + 
Q 



dy) 



p sinh 9 + 



Q 



2! 



f ^ ^ ^ p 2 cosh sinh + 
\dxdyj Q , 



d 2 u 
dy 2 



p 2 sinh 2 9 



Q' 



(59) 



where Q = (q, 0) and Q' is an appropriate point on the segment between Q and the point determined 
byp, 9. 

Let us substitute Eq. (59) in the left hand side of Eq. (57). Some terms do not give contribution, 
because they are anti-symmetric functions integrated in a symmetric range. It results 



+ - 



1 fd 2 u\ 



9x2 )q> 



lim 

29 



p 2 (sinh0*cosh0* + 6*) + 



' [29 i u{Q)}+ 1 



29 l 



1 fd 2 u 



2 v<9r 



p 2 (sinh0 i cosh0 i -0 J 



Q' 



(60) 



Let us express the hyperbolic functions of 9 l as function of p by using Eqs. (55). After this substi- 
tution the limit must be changed as follows lim — > lim . 

As far as the terms in square brackets are concerned in this limit, we have 



• the hyperbolic functions are proportional to 1/p, so the products between powers of p and the 
same power of hyperbolic functions give finite values; 

• the terms p 2 9 l go to zero. 

Then only u (Q) remains and Eq. (57) is obtained. 

The demonstration of Eq. (58) is obtained by means of analogous considerations, taking into 
account that lim p* = 0. But in the final step we divide for the divergent range 2 9 % that is different 

from 2 9 l * at numerator. Now we demonstrate that, in the limit p — > 0, the two divergent ranges 
are equal. 

Actually from Eqs. (24), (26) and (54), let us express cosh#" as function of the same parameters of 
cosh#\ By applying the rule of L'Hospital, we obtain 



lim 



6>*->oo V 29 



29" 



lim 



2 cosh" 1 


p2 


- q 2 + p 2 ' 


2pp 


2 cosh -1 


~p 2 


-q 2 -p 2 ' 


[ 2qp J 



1. 



□ 



(61) 
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Let us extend the demonstration to the general case, in which a / 0. So let us demonstrate, 
referring to Fig. 1, that 



lim 

0' -J-OO 



1 



+6' 



lim 

0« -^oo 



29 l 



29' 

(X*) 



u(p, 9 + a) d9 



Hi) 

u(p*, 9* +a) d9* 



= u{Q) 
= u(Q*) 



(62) 
(63) 



Proof. Let us start from Eq. (62), the Taylor formula of Eq. (59) for hyperbola X, must be 
generalized as follows 



u (p, 9 + a) = u(£ + p cosh (9 + a); rj + p sinh {9 + a)) = u (Q) 



- (Ir) ^c-osh(0 + n)4 

Q 



dy) 



p sinh (9 + a) + 



Q 



2! 



d 2 u\ 



—j /P 2 cosh 2 (9 + a) 



+ 2 



(&a;J e / cosh(9+a)sinh(, ' +a; 



9y 2 



p 2 sinh 2 (6» + a) 



(64) 



where Q = (^, 77) and Q' is an appropriate point on the segment between Q and the point determined 
by p, 9 + a. 

So, in the calculation of definite integrals, sinh 9 l and cosh 9 l are substituted by linear combinations 
of the same functions and all the previous considerations about the limit for p — > hold. 
The same generalization can be done for Eq. (63) and hyperbola X*. □ 



By collecting the results of Eqs. (62), (63) and (49), it results 

u{Q)+u{Q*)=u{P l )+u{P 2 ). (65) 

As a final remark, we give a mathematical meaning to the integrals in Eqs. (57) and (58) divided 
by the diverging angle 29 l . 

Actually, from Eq. (57) we see that the factor 1/(2^) outside the integral, is the same of the 
integration limits, then the left hand side represent the mean value of the function u in the integration 
range. Taking into account Eq. (61), the same meaning can be given to Eq. (58). 
This result is in agreement with the equivalent expressions in the studies of functions of a complex 
variable [8] and for "the initial value problem" for Laplace equation [5]. 

Actually for these problem, for calculating the value of a function in a point P of a domain, given 
its values on the frontier, it is taken a circle around P with radius r — > 0. This circle gives a factor 
l/(2 7r) before the integrals of the right hand side, calculated in the range -H- 2 7r. 



5 Non-omogeneous Wave Equation 

Here we see that the classical approach, by means of integral formulas of Sec. 2.1, allows us to extend 
the obtained results to non-omogeneous wave equation. 

Let v satisfy the omogeneous wave equation and u the non-omogeneous one, that is 

A 2 u = f(x,y). (66) 
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Eq. (9) becomes 

du dv \ 



V dn~ U dn) *~ = J J v f^ x ^) dxd y ( 67 ) 
and, by setting in Eq. (8) v — 1, Eq. (11) becomes 

f^n dr = J J f( x ' y "> dxdy - ( 68 ) 

5.1 Application of Integral Formulas 

Let us do the same steps that are done in Sec. 3.2 for the omogeneous wave equation. Referring to 
Fig. 1, let us apply Eq. (67) to the following domains 

1. domain D between the hyperbolas X and 7; 

2. domain D* between the hyperbolas 7 and X*. 

Let u(x, y) be a function that satisfies the non-omogeneous wave equation (66) in these domains. 
In the application of Eq. (67) to the domain D we set 

v(x, y) — lnr , (69) 

in the application to the domain D* we set 

v(x, y) = v*(x, y) = lnr* , (70) 

where r, r* are given by Eq. (17). 

It can be checked at once that these functions v(x,y) satisfy the wave equation. 
From Eq. (67) we have: 
for domain D 



r ( 



du d\n r\ f Pl ( du (9 In r . 
In r — u — - — ar + / m r — u — - — dr 



'p*( 7 ) \ dn dn J JpU x ) ^ ® n ® n 



J J "(lnr) f(x, y)dxdy , (71) 



for domain D* 



„ du d In r* \ /" p i / 9w 9 In r* \ , 
In r — -u — dr + m r — u — dr = 



P »( 7 ) \ dn dn J Jp^^*) V <9n <9n 



J J (hi r*) f{x, y)dxdy. (72) 



Let us add the left and the right-hand sides of Eqs. (71) and (72) and group together the integrals 
on 7 and the integrals on X and X*. 

From this sum the following total equation results, which includes the integrations on the two domains 
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D and D* 



+ 



hip— dr - 

pi 

+ 



pf 

PI (i) 



In A p — dr 



u 



du 
dn 

din r 

dn 



u 



PI (7) 



d In (r/r*) 
dn 



dr 



dr + 



(r=p) 



/ 

Jpu 



i a du j 
In A p — d r 

dn 



-u 



[ lnr) f(x, y) dx dy + 



u 



PI (I*) 



PI (i*) 

din r 

dn 



1 du A 
lnp dn- dT 



dr 



D* 



(r*=p*) 

In r*) f(x, y) dx dy . 



(73) 



In the equation the left-hand side is the sum of Eqs. (32) and (33), for the same considerations done 
in Sec. 3.2. 

By examining the sum of line integrals in the left-hand side, it results that the third and fifth terms 
make up an integral on the frontier of domain D + D* and the first and sixth terms make up an 
integral on the frontier of domain D*. By applying Eq. (68) to these closed cycles it results 

r-P, 1 a.. rPi 



f 1 du ft du f f 

/ lnp—dr+ / lnp — dr = lnp / f(x,y)dxdy, 

Jpid) an J Pi fx*) an J d+d* J 

f(x, y) dxdy 



I In Ap — dr 

J PI ( 7 ) on 



pi 



1 A d U 7 

In A p — — d r 

dn 



In Ar. 



D* 



(74) 
(75) 



and, Eq. (73) becomes 



r - 

J PI (X) 



din r 

dn 



[Pi 

dr + J a 

(r= P ) J Pi (x*) 



d In r* 

dn 



dr 



(r*=p*) 



U 



PI (7) 



d In (r/r*) 



dr 



= — lnp 



/ / /(*, 

J D+D* J 



y) dx dy + In A p 



D* 



dn 

f(x, y) dxdy 



+ J j "(lnr) f(x, y) dx dy + J J ( In r*) f(x, y) dx dy . 
The last equation is the extension of Eq. (34) to the case of non-omogeneous wave equation. 



(76) 



5.2 Limits of Area Integrals 

In Sees. 4.2 and 4.3 the terms of the left-hand side of Eq. (76) are divided by the integration range 
2 6 l where, from Eq. (55), 

(77) 



9 l = cosh 



-1 



2qp 



and the limit of the ratio is calculated for p — > (i.e., 1Q % — > 00). In this way the results of Eqs. 
(62), (63) and (49) are obtained. 

Here the same procedure is applied to the terms of the right-hand side of Eq. (76), in order to extend 
the Eq. (65) to the case of non-omogeneous wave equation. When p — > the hyperbolas X and X* 
become the parallels to axes bisectors from Q and Q* (see Sec. 2.2 and Fig. 1) and the domains D 
and D* are contained by the arc P1P2 of hyperbola 7 and the straight line segments QPi, QP2 and 
Q*Pi, Q*P2, respectively. 

The following limit values result for the four area integral terms of Eq. (76) 
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1. 



lim 



-inp / / f(x, y) 

J D+D* J 



dx dy 



2 cosh" 1 



p 2 — q 2 — p 2 



2qp 



D+D* 



f(x, y) dxdy. 



(78) 



Proof - The result of Eq. (78) is obtained by applying the rule of L'Hospital to calculate 



lim 

p^-0 



In p 



2 cosh" 1 



p 2 - q 2 - p 21 



2qp 



1 

2 



(79) 



Also it may be obtained by substituting the expression 



cosh 1 



p 2 — q 2 — p 2 



2qp 



= In 



p 2 — q 2 — p 2 



2qp 



+ 



1 1 2 \ 2 

p z — q z — p zx 



2qp 



- 1 



(80) 



and calculating the limit directly. 



□ 



2. 



In A r 



lim 



D" 



f(x, y) dxdy 



2 cosh" 1 



p 2 — q 2 — p 2 



2qp 



Proof - Eq. (81) is obtained by the direct limit calculation 



lim 



In A r 



2 cosh" 1 



p 2 - q 2 - p 2 



0. 



3. 



lim 

p^O 



2qp 

J J O^r) f(x, y) dxdy 



2 cosh" 1 



2 2 2 

p z — q z — p A 



2qp 



(81) 



□ 



(82) 



0. 



(83) 



Proof - Let M be the absolute maximum of \ f (x, y)\ in the domain D at limit p — > 0. Let us 
transform the coordinates x, y into polar coordinates 



x = £ + r cosh(# + a) ; 
The area element is transformed 



y = i] + r sinh(# + a) . 



It results 




dxdy = r dr d9 . 



lnr) f(x, y) dx dy 



ID 



< M 



(\nr)rdr d6 



(84) 
(85) 

(86) 
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where ±6 r are the extreme values of 6, corresponding to the intersections points of the hy- 
perbola that has center in Q and semi-diameter r, with hyperbola 7. The hyperbolic angled r 
is given by Eq. (77) or (80), where p is substituted by r. Then, by using the logaritmic form 
(80), 




lnr) f(x, y) dxdy 



'D 



2M 



(p-q) 



lnr) r < In 



< 



p 2 — q 2 — r 2 



2qr 



+ 



2 2 2 \ 2 
_ _ r ^ \ 



2qr 



dr 



(87) 



The function in the integral in the right hand side has finite values for p < r < (p — q) and, 
when r = p and p — y 0, 



lim < ( In r) r In 



r-s-0 



p 2 — q 2 — r 2 



2qr 



+ 



2 2 2 \ 2 
_ _ r ^ \ 



2gr 



— >■ lim < ( lnr) r In 



r-i-0 



p — g 
qr 



lim {( lnr) r (- lnr)} = . 



(88) 



The result of Eq. (88) is obtained by applying two times the rule of L'Hospital. Therefore the 
function in the area integral at numerator of (83) is finite in all the domain D, also in the limit 
p — > 0. Then the corresponding area integral is finite and, divided by the hyperbolic angular 
range that diverges when p — y 0, gives the result of Eq. (83). □ 



4. 




In r*) f(x, y) dx dy 



2 cosh -1 



p 2 — q 2 — p 2 



2qp 



0. 



(89) 



Proof - A process of demonstration analogous to the one developed for Eq. (83) can be used, 
recalling that lim p* = . □ 

Finally, by collecting the results of Eqs. (62), (63), (49) and (78), (81), (83), (89), the following 
equation is obtained 



u(Q)+u(Q*)=u(P 1 )+u(P 2 ) + \ [ [ f(x,y)dxdy, 

1 Jd+d* J 



(90) 



that is the extension of Eq. (65) to the case of non-omogeneous wave equation. 



6 Conclusions 

By studying the wave equation in a plane with its own geometry, similar results to the ones of Laplace 
equation studied in Euclidean plane are obtained. 

These results can also be read in the following way: in [1] the Euclidean theorems have been used 
as starting points for their extension to hyperbolic geometry by means of analytical demonstrations. 
In this paper, this method is extended to a problem, related to wave equation that is usually studied 
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in Euclidean geometry. 

We know that the wave equation is the starting point for obtaining the Lorentz transformation 
[3] from which a physical meaning to hyperbolic geometry is given. Therefore the application of 
hyperbolic geometry can be considered a natural way to study the wave equation in a Cartesian 
plane. 

Taking into account the theoretical and practical relevance of wave equation in Mathematics and 
Physics and the many subjects related to the treated arguments, it can be expected that the obtained 
results may be the starting point for improvements in many directions. 



A Normal and Tangent Local Coordinates in the Hyper- 
bolic Plane 

Here we construct the local reference frame used for calculating line integrals in a hyperbolic plane. 
This reference frame has the origin in the points that move along the curve and an axis tangent to the 
curve. It is equivalent to the local reference frame employed in the line integrals in Euclidean plane. 
Otherwise the different topology and metric properties [1] generate different geometric characteristics. 
Calling C(x c ,y c ) the generic point on the curve, a local Cartesian reference frame is defined with 
origin in C. The coordinates on the local axes, parallel to x and y axes, are x — x c and y — y c . 
In this reference, let us consider another reference, taking a r axis tangent to the curve and oriented 
according with the integration direction. It forms a hyperbolic angle a with respect to y local axis. 
The other axis, n, forms the same hyperbolic angle a with respect to x local axis, n and r axes are 
orthogonal in the hyperbolic geometry, i.e., they are symmetric with respect to the axes bisectors of 
the local Cartesian reference frame. These axes are oriented so that the n, r frame is congruent with 
the x, y frame. 

In Fig. 3 an equilateral hyperbola is represented and four positions of the point C are considered. 
For each position the two possible directions of r axis are reported. In this way the types of pairs of 
unity vectors defining n, r reference axes are described. 

The transformation equations that link the coordinates x, y to the coordinates n, r are 

x — x c = ± [n cosh a + r sinh a] 

y — y c = ± [n sinh a + r cosh a] (91) 

with the inverse one 

n = ± [(x — x c ) cosh a — (y — y c ) sinh a] 

t = ± [— (x — x c ) sinh a + (y — y c ) cosh a] . (92) 

In these equations the sign + is used when r is oriented upward and the sign — when r is oriented 
downward. 

Let us consider a function u defined in the plane x, y. By using Eqs. (91) and (92) let us express, 
as a function of the local variables n, r, the differential form 

du , du , . . 

3i iv+ 3i ix (93) 

By considering a line integral of this form, it results dn = 0, so that 

dx dx dx 
ax = — dn + — — dr = -^-dr 
on or or 

dy = tt dn + ^ dr = ^ dr . (94) 
on or or 
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Figure 3: Normal (n) and tangent (r) unity vectors along a hyperbola in the hyperbolic 
plane. 

The topology and geometry of hyperbolic plane generate different definitions from the ones of Eu- 
clidean geometry, now we recall [1] and represent the orthogonal lines. 

We call a the angle between the r axis and the parallel to y axis of the local reference frame. 

In this figure we show how normal (n) and tangent (r) unity vector pairs are defined, with respect 

to t orientation, in the points of both the arms of a hyperbola and for y > or y < 0. 

The solid lines represent r unity vector in the upward direction. The n unity vector is in the right 

sector of the local reference frame, a is positive when r unity vector is on the right hand side with 

respect to the positive local Cartesian y axis and is negative when it is on the left hand side. 

The dashed lines represent r unity vector in the downward direction. The n unity vector is in the 

left sector of the local reference frame, a is positive when r unity vector is on the left hand side with 

respect to the negative local Cartesian y axis and is negative when it is on the right hand side. 
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On the other hand, from Eq. (91), the following relations between the derivatives with respect to 
the orthogonal directions of the tangent and the "hyperbolic normal" to the curve result 

Ox Oy Ox Oy 

— — = ± sinh a = — — — — = ± cosh a = — — (95) 
Or On On Or 



and the differential form of Eq. (93) becomes 

Ou Ou f Ou Oy Ou 0x\ _ (Ou Ox Ou 0y\ Ou . . 

Ox ^ Oy X \0x Or Oy Or J T \0x On Oy On J T On ° 



Partial Derivative with Respect to Normal Local Coordinates. 

Referring to Fig. 3, partial derivatives of x and y with respect to n are calculated with the 
conditions 

• the point of coordinates x, y is on the equilateral hyperbola 7 with center in O and semi- 
diameter p. Its equation can be represented by Eq. (14), by setting the line parameter <fr+a = a. 

• t is oriented upward. 

From Eqs. (91) in which + sign is used and Eqs. (14), we have 

Ox x Oy . , y , n _ s 

— — = cosh a = — ; — = smh a = -. (97) 
On p On p 
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